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The quantum mechanical counterpart of the famous Stoner-Wohlfarth model - an easy-axis mag¬ 
net in a tilted magnetic field - is studied theoretically and through simulations, as a function of 
the spin-size S' in a sweeping longitudinal field. Beyond the classical Stoner-Wohlfarth transition, 
the sweeping field-induced adiabatic change of states slows down as S increases, leading to a dy¬ 
namical quantum phase transition. This result is described as a critical phenomenon associated 
with Landau-Zener tunneling gaps at metastable quasi-avoided crossings. Furthermore, a beating 
of the magnetization is discovered after the Stoner-Wohlfarth transition. The period of the beating, 
obtained analytically, arises from a new type of quantum phase factor. 

PACS numbers: 75.10.Jm, 75.45.+j, 75.50.Xx, 75.60.Jk 


We study the reversal of a uniaxial magnet of spin 
S submitted to a fixed transverse field and a longitu¬ 
dinal sweeping field. This model had been studied for 
single-molecule magnets with moderate spins 5'=10, in 
which the Landau-Zener (LZ) transition plays an impor¬ 
tant role [iHi- Some years ago, the parent model of 
an Ising spin-chain with ferromagnetic interactions and 
same longitudinal and transverse fields showed a quan¬ 
tum spinodal phase transition where a size (spin-chain 
length)-independent magnetization decay was observed 
and attributed to independent single-spin reversals 

In this letter, we study the quantum aspects of the 
Stoner-Wohlfarth (SW) transition, and more particularly 
the relations between the quantum dynamics and classi¬ 
cal irreversibility of a large quantum spin S' in a uniaxial 
anisotropy broken by a transverse field, and submitted 
to a sweeping longitudinal field at zero Kelvin. In the 
S ^ oo limit, the spin becomes classical and exhibits the 
so-called SW transition Q. This is a magnetization jump 
from a metastable to a stable state when the field, ap¬ 
plied in the opposite hemisphere, reaches a critical value. 
We obtained two main results: (i) the classical SW tran¬ 
sition (infinite spin) is given by a critical phenomenon of 
the spinodal type in the limit of a large spin S in the 
quantum regime, and (ii) when the sweeping field ex¬ 
ceeds the SW point, magnetization beatings of quantum 
mechanical origin, which are analyzed and attributed to 
a new type of quantum phase factor, appear. 

In order to catch the properties of our model prop¬ 
erly in the S ^ oo limit, we introduce the normalized 
quantum spin operators with a modified commutation 
relations: 


S i 

{a = X,y,z), [Sa,Si3] = -eal^-yS-y. (1) 


The corresponding SW Hamiltonian, with uniaxial 
anisotropy, transverse field (fixed) and longitudinal field 


(sweeping at the time-rate c), is written as 

n = -Dsl - - H,s,, H, = iff) - ct. (2) 

In Eq. (|2]) and hereafter, we set gfiB = I and h= 1. 

The time evolution of the normalized quantum spin 
operators, given by 

dsjQ I ^ ^ 

~ ^ \^y^Z + SzSy) + HzSy} , 

ds I 

{ D (Sz^x T ^x^z) HzSx T HxSz\ , (3) 

dsz _ I ^ 

is obtained numerically by the standard Runge-Kutta 
method. 

The corresponding time-evolution of the usual (classi¬ 
cal) SW model comes from the torque equation dm/dt = 
—mx iTeff with the effective field iTeff = —dEsw/dm = 
{Hx, 0, 2Dmz -P Hz). The energy of the SW model, Esw 
, is of course derived from Eq. 0. Comparing with Eq. 
m, we find that the dynamics becomes the same if we 
ignore the commutation relations among spin operators 
in Eq.(j3]). But the time in the quantum system should 
be normalized as 


T = t/S, (4) 

to study correspondence to the classical dynamics. 

Motion of a spin of classical SW model m(t) drows a 
trajectory on the unit sphere. If the fields Hz and Hx are 
weak, due to the anisotropy D, the system has metastable 
state in which the Sz is antiparallel to the direction of 
Hz, but the state is locally stable. If Hz and Hx increase 
and the SW condition (211)^/^ = (iL^,)^/^ + {Hz)‘^^^ is 
satisfied, the metastable state becomes unstable. 

Fig. □ shows the classical motion of magnetization 
{mx{t),my{t),mz{t)) under a sweeping longitudinal field. 
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FIG. 1: (color online) A classical trajectory under the sweep¬ 
ing field Hz = 4: — ct {c = 0.01) for D = 1, = 1. Arrow 1 

indicates the starting point which is a metastable fixed point 
for Hz = 4, and the arrow 2 the point of the SW jump fol¬ 
lowed by precession around the magnetic field. 



FIG. 2: (color online) Energy spectrum for S' = 20 with D = 1 
and Hx = 1. as a function of Hz. For each value of Hz 
there are 25' + 1 = 41 eigenvalues which are plotted by blue 
points. The brown line shows the position of the SW points. 
The red bar denotes the time evolution of population at each 
eigenstate under a sweeping field Hz = I — O.OSt. 


It starts from the metastable point at Hz{0) = 4 (arrow 
1 at t = 0) and ends at the SW point (arrow 2 t = tsw)- 
The motion in this process is adiabatic. At t = tsw 
the irreversible magnetization jump takes place. It is 
followed by a precession-like motion about the local neg¬ 
ative effective field, the amplitude of which decreases as 
Hz{t) becomes larger. 

To characterize the SW transition of a large quantum 
spin 5, we first investigate the energy-level structure as a 
function of Hz. The example 5 = 20 given in Fig.[2]shows 
how the 25+1 levels with positive and negative slopes 
(i.e. spins), intercept at avoided-level crossings with gaps 
/AE. If a sweeping field is applied to the ground-state 
(positive side, lowest line), say from Hz = 1, this state 
remains occupied until Hz = 0 where the ±5 avoided- 
level crossing is reached. If the reciprocal sweeping-time 
is much smaller than AE at this point, the motion is 
almost adiabatic and the spin-reversal probability is close 


to unity (LZ model [l|). However, in practice this change 
hardly occurs because, with realistic parameters, AE is 
vanishingly small for a large 5 and only a very small 
fraction of the population scatters from +5 to —5. The 
rest remains on the line of slope 5, continuing the Hz > 
0 ground-state into the Hz < 0 region. This Hz < 0 
line corresponds to the classical metastable state and can 
be called the metastable braneh (See also Fig. [2j). This 
branch crosses the levels of slopes M = —5 + 1,—5 + 
2,''' ,—5 + /c,''' successively. At each crossing, some 
population of the metastable branch scatters to the state 
of negative magnetization (Fig. [2]). The Landau-Zener 
probability of the population remaining at the metastable 
state, after the k-th avoided level-crossing is given by 


(S) 

Pk = exp 


^ 7r(A+))A 

4fic(®)AM+ 1 ’ 


( 5 ) 


where is the energy gap at the field and, result- 

ing from the mixing of the spins 5 (metastable branch) 
and —5 + /c (intercepted branches), AM^^^ ^ {2S — k)/S 
is the difference of magnetization (slopes) of the two 
states. As 5 becomes larger, the spectra associated with 
a normalized spin densify because the number of the 
eigenvalues (25+1) increases. For example, going from 5 
to 25 leads to a new eigenvalue between two consecutive 
eigenvalues of the initial spectrum of 5. In the contin¬ 
uous limit, the amount of scattering within a given-field 
interval must be nearly same, leading to the condition 
P^ 2 k^P^ 2 k+i ~ P^k^ successive LZ transition proba¬ 
bilities. Due to the spin renormalizations {sz = Sz/S, 
Eq. ([I])), one must have = AMjf^ and ct = vr 

with V = c5, so that the realization of the probability 
condition given above implies rescaling the gap, and AE 
becomes SAE. The latter is plotted vs. Hz (Fig. [3|). 

Finite AE means that the magnetization can reverse 
at fields \Hz\ < |i^sw| with some probabilities, whereas 
in the 5 ^ oo limit, the gap is null and the reversal 
is prohibited as expected in the classical limit. This is 
simply a quantum tunneling effect in the SW model. For 
example, the quantum SW transition begins to take place 
at Hz ^ 2iLsw/3 for 5 = 20 (Fig. [3]). The dependence of 
the scaled gap SAE on the spin changes at Hz = 
and it remains finite at \Hz\ > |iLsw|: which indicates 
that the spin undergos scattering to state of positive mag¬ 
netization. This apparently surprising result is simply 
due to the non-adiabatic character of the SW transition. 
The similarity of the field variations of the scaled gaps of 
Fig. [3] suggests some kind of criticality which will now be 
rapidly investigated. 

In the mean-field approach, the spinodal phase transi¬ 
tion of a classical system of size N, follows the scaling plot 
r = N^^^f{{H - where r is the relaxation 

time near the spinodal point H = Hsp 0. If our total 
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FIG. 3: (color online) Normalized energy gaps at avoided- 
level crossings plotted vs the sweeping magnetic field for S = 
20 (green), 40 (blue), 80 (purple), 160 (cyan), and 320 (red) 
with Hx = 1, D = 1. (inset) Scaling plot showing a mapping 
of the LZ transition to a dynamical spinodal phase-transition. 





FIG. 4: (color online) Beating observed after the SW point. 
Note that the data are plotted as a function of Hz because Hz 
is swept linearly in time i.e., Hz{t) — Hz(fi) — ct. The red line 
shows the fidelity and each green, blue, and orange 

lines show (s^), {sx) — 1-5, and {sy) — 2.5 where parameters 
are S — 20, Hx = 1, D = 1, and v — 0.08. 


spin S could be regarded as the size of a spinodal sys¬ 
tem, then the data of Fig. [3] should be plotted according 
to the scaling form, r = S^^^f{{H — Hsp)S‘^^^)^ where 
r oc exp(7r(5'F^^‘^^)^/4/ic^‘^^In Fig. [H inset), we 
find that the data collapse well in a same curve, giving 
clear evidence of the spinodal character of the quantum 
SW model. 

We now switch to the last part of this letter, relative 
to magnetization beatings. When the applied field is 
swept from right to left in Fig. O the metastable curve 
is reached in negative fields inducing population scat¬ 
terings, across avoided-level crossings, taking place from 
the SW point. It is in this region that the magnetization 
beatings are observed (Fig.|4]). We define the spin-length 
fidelity as s(t)^ = {sx{t))‘^ + This 

quantity is not conserved, in contrast to the spin modu¬ 
lus 5'(5'+ 1). The time evolution of each spin-component 
and of fidelity is plotted in Fig. IH The fast oscillations 
are simple precessions (as in Fig.[T]) whereas the beatings 
are modulations of these precessions. The three magneti¬ 


zation components show the same beating period, which 
means that the spin-length s{t) is also beating with this 
period. The amplitude of the Sz{t) beatings decreases as 
the field becomes large, because large fields induce a tilt 
of the precession plane towards the x^-plane, while the 
xy components keep the same amplitude. Contrary to 
precessions, beatings are not seen in the classical case. 
They can be viewed as a modulation of the precessions 
of genuine quantum origin. As this will be shown now, 
they originate from a new type of quantum phase factor. 
The time evolution of the sweeping-field state is given by 
the time-dependent Schrodinger equation id^ jdt = 
where H is the scaled Hamiltonian (j2j) with H^z^ = 0 for 
simplicity. Introducing the unitary transformation to a 
non-linear acceleration rotation frame 


— with U = exp 

it comes 

.d 


i ( Dslt - ) , (6) 


= -HxUUxU4> = 


with 


0(r) = 2 '^^^ ~ 2Dszr. 


( 7 ) 


Urot =-Hx{sx cosQ{t) - SySmQ{T)}e ( 8 ) 

where 


( 9 ) 


It consists of a Zeeman term in the chosen rotating 
frame with classical phase given by expression mul¬ 
tiplied by a quantum phase term (or in 

S units). Note that the sweeping velocity v = cS = ctjT 
intervening in Eq. (jS]) results from Eq. 0. The first term, 
which is independent of S', tends to the classical model in 
the S ^ oc limit. It represents the classical fast preces¬ 
sion motion of Eig. [H The second term depends 

on S and tends to unity in the S ^ oc limit. It is of gen¬ 
uine quantum origin and is responsible for the observed 
beatings. Its specific field-period is given by 


Th^ = cTt = c 


27r 

D/S^ 


2ttvS 
D ’ 


( 10 ) 


where Tt is its specific time-period. The period is propor¬ 
tional to the spin size S', the classical sweeping velocity 
and inverse proportional to the anisotropy constant as 
observed in the simulations (Eig. [5]). This confirms that 
the ’’beating” phenomenon is due to the factor . 

Einally, we should note that if the field sweeps back, this 
quantum beating disappears of course when the field goes 
through the strongly irreversible SW point. Interestingly, 
in the reversible region above the SW point, the beatings 
persist, even if the field is stopped or cancelled at a given 
time. All that is illustrated in Eig.[6]after the third beat¬ 
ing, when the field is kept constant (Eig. [6] (left)) or is 
increased again (Eig. [6] (right)). 
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FIG. 5: (color online) Dependence of the ^-component of 
magnetization dynamics for different values of: (left) spin S = 
20 (red), 30 (green), and 40 (blue) for v = 0.08 and D = 1.0, 
and (right) anisotropy constant D = 1 (red), 2 (green), and 
3 (blue) for S' = 20 and v = 0.08. The dependence on v 
(v = 0.08, 0.12, and 0.16 gives the same figure as (left). 



FIG. 6: (color online) Beating of the spin-length fidelity, and 
of the three spin components in a sweeping field above the SW 
transition for S = 20 for D = l^Hx = 1 with the sweeping 
velocity c = 0.004 until Hz{t) = —12. (Left) Beating is not 
altered after the field (brown curve) is set constant. (Right) 
Beating is suppressed after the field (brown curve) is increased 
again and re-enters above the |iLsw|- In these figures, we plot 
the time evolution as functions of time but not Hz{t) as in 
the previous figures. 


In conclusion, we have studied the dynamics of the 
classical and the quantum SW models, showing sur¬ 
mountable and impassable bridges between the two and 
their identification when S ^ oo. The two most impor¬ 
tant results consist in (i) the description of the quantum 
SW model in terms of a dynamical spinodal phase transi¬ 
tion with a scaling of the tunneling gap (or, equivalently, 
of the time associated with the tunneling probability) vs. 
Hz — Hsw- A detailed description of a population scat¬ 
tering effect (to the states of positive or negative magne¬ 
tizations) occurring across a succession of LZ transitions 
along a continuation of the positive-field ground-state 
in the negative fields above the SW transition (that we 
called metastahle branch). When the field sweeps below 
the SW transition, the gaps - and therefore the spin- 
reversal probabilities - at avoided-level crossings are so 
small that the system almost remains metastable, but 
dynamical transitions to the stable point is still possi¬ 
ble. However, when S ^ oo the dynamics disappears 


leading to the classical, static, SW model at T = 0. A 
spinodal scaling (Fig. [3j) gives a synthetic representation 
of this complex physics, paving the way for more general 
studies on the quantum to classical transition. 

This rich physics associated with the metastable 
branch also leads to our second important results (ii) 
when the sweeping field is larger than the critical SW 
field, the spin motion, classically described by a time- 
dependent precession about its slowly moving local-field, 
becomes modulated in time leading to a beating of the 
three spin components and of the spin-length fidelity with 
the characteristic period 2'kvSID obtained analytically 
thanks to a unitary transformation allowing one to de¬ 
compose the classical motion and the quantum mechan¬ 
ical phase terms, by which the period of the beating is 
well explained. 
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